Current- Voltage Characteristics in Donor- Acceptor Systems : Implications of a 

Spatially varying Electric Field 
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We have studied the transport properties of a molecular device composed of donor and acceptor 
moieties between two electrodes on either side. The device is considered to be one-dimensional 
with different on-site energies and the non-equilibrium properties are calculated using Landauer's 
formalism. The current-voltage characteristics is found to be asymmetric with a sharp Negative 
Differential Resistance at a critical bias on one side and very small current on the other side. The 
NDR arises primarily due to the bias driven electronic structure change from one kind of insulating 
phase to another through a highly delocalized conducting phase. Our model can be considered to be 
the simplest to explain the experimental current-voltage characteristics observed in many molecular 
devices. 



The study of electron transport through single 
molecules is gaining tremendous attention in recent years 
owing to the wide variety of applications that they can 
be used in^ . Recent advances in experimental techniques 
have allowed fabrication and measurement of current 
through such nanoscale systems. Various molecular sys- 
tems have already been demonstrated to behave as wires, 
switches, diodes and RAMsSiSii. The ability of a molecule 
to switch between off and on states is one of its most im- 
portant applications. Experimentally, this has been ob- 
served in several organic molecules with various donor 
and acceptor substituents^iSi^. Considerable amount 
of theoretical work based on semi-empirical to a6-initio 
methods have also been performed to model molecular 
transport characteristicsi. Many explanations for this 
switching phenomenon, based on charging2iifl, reduction 
of the acceptor moietji^, twisting of the ring structure 
leading to conformational changesiiiiSiiii^, bias driven 
changes in molecule-electrode couplingi^, have been pro- 
posed. Most of these require to impose some external 
factor like the rotation of the middle ring etc in order 
for the external bias to cause NDR at some bias. And 
most often they do not make a relation between struc- 
tural preference and bias polarity and hence do not ex- 
plain the asymmetry that has been observed in the ex- 
periments. Also a comprehensive understanding of the 
switching phenomena in general cases is still elusive. 

In this letter, we try to understand the reasons for 
the observed Negative Differential Resistance (NDR) in 
molecular wires based on a very simple donor-acceptor 
model. In fact, the asymmetry as well as NDR in for- 
ward bias that has been observed in the Tour moleculesS*^ 
comes out naturally out of this dimer model without hav- 
ing to invoke any external factors. And the simplicity 
of this parametrized model makes it very tractable and 
lends physical insight into the factors causing the NDR. 
We find that a spatially varying external forward bias 
switches the electronic phases resulting in a sharp rise 
and fall in transmission through the device. Various ex- 
tensions of the dimer model together with different spa- 
tial variations of the electrostatic potentials have also 



been considered. 

A closer look at the structure of the Tour molecules 
(see inset of Fig.l) suggests that one part of the molecule 
with the donor group (NII2) has a positive on-site energy 
and the other with acceptor group (NO2) has a negative 
on-site energy (Although this has been assumed here, we 
have later extended it to describe systems with many lo- 
calized segments in it, as is shown by quantum chemical 
calculations on the Tour molecules). For such a two-level 
model involving a donor and an acceptor, the Hamil- 
tonian is given as: H = 2^i'A^i + ^(^102 + he), 

where is the on-site energy of site i and t is the hop- 
ping integral. This dimer is attached to two metal- 
lic electrodes on either side, which are assumed to be 
non-interacting semi-infinite one dimensional systems de- 
scribed by a tight-binding Hamiltonian with an energy 
bandwidth of 47 where 7 = lOeV—. The Fermi energy 
of the electrodes is generally assumed to be a fitting pa- 
rameter. We assume that the Fermi energy lies halfway 
between these equilibrium (zero bias) energies, though 
its location does not influence our main results. 

Spatial variation of bias on the structure is quite a 
central issue in this fieldi^. Except for a few cases (given 
later), we consider that it drops on the device as a ramp 
function, varying linearly from one electrode to the other 
as: Vn — —nV/{N + 1), where n is the site-index and 
N, the total number of sites. With the potential, the 
energies for the dimer can be written as: 
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The coupling to the electrodes modifies the bare Greens 
function of the molecule, which across the molecule can 
be written as 
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FIG. 1: The current-voltage characteristics for the 2-site sys- 
tem for £2 = — £i = 0-5eV and t = O.leV^. Inset is the Tour 
molecule. 

where Ei and S2 are the self-energies corresponding 
to the two electrodes, calculated within the Newns- 
Anderson modelAS-. Using the Greens function, the cur- 
rent through the system can be obtained from the Lan- 
dauer's formula^S: 

HV)^^r' dE[Tr{TiGT2G^)] (3) 

JEf-eV 

where ri^2 are the anti-hermitian parts of the self-energy 

matrices, ri_2 = *(5^i,2 — 5]| 2) which describe the broad- 
ening of the energy levels due to the coupling to the elec- 
trodes. 

In zero bias condition, as can be seen from Eq.l (with 
V = OV), presence of different on-site energies opens 
up a gap larger than that for a purely hopping model 
(£2 = El = OeV) near the zero of energy indicating the 
preference of the electrons to stay at the atomic site with 
negative on-site energy. The equilibrium transmission is 
found to be large for purely hopping model since it cor- 
responds to equal distribution of charges. With the in- 
clusion of different on-site energies, the systems becomes 
insulating due to charge transfer and the zero-bias trans- 
mission reduces due to this preferential charge localiza- 
tion. 

Fig 1 shows the nature of the current-voltage charac- 
teristics with external bias. As can be seen, the current is 
negligible around the zero of energy and around a bias of 
IV, there is a small jump in the current. This is the bias 
at which the device energy level comes into resonance 
with the Fermi energy of the electrodes-^. With increase 
in the forward bias, around a bias of 3^, the current 
shows a sharp rise and fall, indicating strong Negative 
Differential Resistance (NDR). On the other hand, with 
increase in the reverse bias, the system continues to re- 
main insulating with negligible current. 

To understand the reasons for the NDR, we look at the 
variation of the energy levels (Ek) of the bare molecular 
dimer with bias (Fig 2a) and the numerator of the Greens 
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FIG. 2: (a) The variation of the two levels (circles and stars) 
of the 2-site system with the applied bias, for £2 ~ — ei = 
0.5eV, t = O.leV The dotted lines indicate the variation of the 
Fermi energies of the electrodes with bias, (b) The numerator 
of the Greens function matrix element for the corresponding 
energy levels shown in (a). 1 and 2 represent the site index 
and k specifies the corresponding level, (c) The IPRav for 
the levels shown in (a). See the text for definition. 

function, (l|fc)(2|fc) (see Fig 2b, fc = 1,2 are the eigen- 
states). With increase in the forward bias, the energy 
levels come close to one another up to the critical bias Vc 
at which the NDR is seen, above which, they move far- 
ther away. In Fig. 2b, exactly around this Vc, the contri- 
bution to the eigenstate (MO) coefficients from the sites 
increases quite sharply. To quantify this critical bias Vc, 
we minimize the gap between the energies with respect 
to the applied bias and obtain, Vc = 3(e2 — ei), which is 
in accordance with our numerical data (for ei = —O.beV 
and £2 = 0-5eV, we find Vc ~ 3^)^^. At this critical bias, 
the energies take the values =Fi, precisely the energies of 
the non-interacting system. However, with increase in 
the reverse bias, the energy levels start diverging away 
from their zero bias gap making the system more and 
more insulating, explaining the small current that is ob- 
served in Fig.l for negative bias. 

To quantify the effect at Vc, we calculate the Average 
Inverse Participation Ratio {IPRav) which defines the 
extent of localization for a given eigenstate, with Ek'. 

where P'^ is the IPR, defined as P'^ = ^ J2j IV'(j, k)\^ 
where the j is the atomic site index and D{E) is the 
density of states. Fig. 2c shows a strong dip in the values 
of IPRav around the critical bias indicating complete 
delocalization in the system, while at other values of the 
bias, IPRav is much larger due to the localized nature of 
the eigenstates. 

Initially for small bias, as noted before, the system 
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FIG. 3: The I-V characteristics of a 20 sites system, for £2 = 
— ei = 0.5eV and t = l.OeV. The inset shows the variation 
of the number density at the even sub-lattice for the case 

i (triangles) and case ii (diamonds) discussed in the text. 
Similarly in the odd sub-lattice for case i (circles) and case 

ii (stars). Also shown the pictorial representations of the 
corresponding eigenfunctions in site basis, very close to the 
critical bias. Filled circles indicate large contribution to the 
eigenstate. 



tends to accumulate its charge density at the site with 
lower on-site energy. Such a localization causes the sys- 
tem to be insulating. If this site is closer to the electrode 
with higher chemical potential, charges tend to move to- 
wards the other site. At Vc, where the NDR is seen, 
the charge densities are equally distributed at both sites 
with no preference of one site over another indicating a 
situation where both the on-site energies are equal. Fur- 
ther increase of bias would localize the charges on the 
other site resembling an insulating dimer with its on-site 
energies interchanged, precisely the case as with reverse 
bias. 

To extend these dimer results further, we look at a 
chain of N atoms with alternating on-site energies, ei 
and 62 . The energies in the absence of external bias have 
the form: 



E. 



ei -e2)2 + 16t2cos2( 



N + 1 



)] (5) 



where N is the number of sites and the '—(+)' signs 

hold for j = 1 ,2 . . . . , -J (j = Y + 1 ' N) . It is quite well-known 

that different diagonal terms (either random or alternat- 
ing) localizes the eigenstates of a one-dimensional sys- 
tem. From a perturbative treatment, this localization 
can be made quantitative as 1/A = log where 
the localization length A depends on the width of the 
diagonal term W and the hopping t. Since the sys- 
tem becomes highly delocalized at we can naively 
expect the width to tend to zero at this bias, giving 
Vc\t^o = (£2 - ei){N + 1)/{N - 1), for a N sites chain. 



For a 20 sites half-filled system, we obtain a K of \.GV 
for values of ei — — O.Se^ and £2 = Q.beV and t = leV, 
which is close to the value calculated from the above ex- 
pression (~ I-IV), with second-order corrections due to 
the hopping termSi. Fig 3 shows the I-V characteris- 
tics for this system which is asymmetric and shows clear 
eigenvalue jumps for reverse bias^^ and a sharp NDR 
peak at a critical forward bias. 

To quantify the effect of bias in the system, we cal- 
culate the total charge densities in each of the two sub- 
lattices for two cases: (i) charges are filled up to the 
highest-occupied level (HOMO) and (ii) with one of the 
highest occupied level electron promoted to the lowest 
unoccupied level (LUMO). Both of these, for each sub- 
lattice are plotted in the inset of Fig. 3, as a function of 
external bias. As can be seen, there is a jump in the 
charge density in the +e sub-lattice due to the initia- 
tion of charge transfer from — e sub-lattice to the sub- 
lattice with +e, exactly at the critical bias. Since the 
HOMO and LUMO levels are very close in energy at the 
critical bias (as shown for the 2-level system in Fig. 2), 
the charge densities corresponding to both cases become 
equal for each sub-lattices, giving rise to the crossing of 
the respective curves. We have also shown in the figure, a 
pictorial representation of the HOMO and LUMO eigen- 
states close to Vc- As can be seen, the contribution to 
HOMO primarily comes from one end of the chain and 
to LUMO from the other end. However, since at Vc these 
two states are quasi-degenerate, a linear combination of 
the two correspond to a system with large transmission 
amplitude between two ends. 
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4: The I-V characteristics for the 20 site system with 
— ei = 0.5eV and t — leV with ramp potential (dia- 
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monds), 1=2 drop(plus), and 1=8 drop (circles) close to the 
interface. The inset shows the variation of the HOMO-LUMO 
gap with bias for the three cases. 

As is well-known, due to screening effects, the potential 
gradient is expected to be larger near the electrodes than 
at the center of the chain where it may be zero or very 
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smalP^. We have considered the spatial variation of the 
potential over / sites close to the interface with variation 
in I, and the potential is a ramp function when I = N. 
Three different potential drops have been considered: (a) 
ramp (h) I — 2 (2 sites drop) and (c) Z = 8 ( 8 sites drop) 
close to the interface. In Fig 4., we have plotted the I- 
V characteristics for the system with £2 = — ei = 0.5eV 
and t = l.OeV for all the three cases. For the one with 
I — 2, the I-V curves shows clear eigenvalue jumps. With 
the increase in I, the NDR begins to appear at smaller 
voltages, and a very sharp NDR peak is seen for the ramp 
potential. The inset shows the variation of the HOMO- 
LUMO gap with bias for the three cases. There is no 
significant closing in of the levels when 1 = 2. But with 
increasing I, this becomes more and more pronounced 
and the NDR appears earlier and is sharper. This study 
indicates the importance of the spatial variation of the 
applied electric field in order to see an NDR response. 

In the long-chain limit, there is an implicit relation 
between variation in on-site energies and the bond-length 
alternation (BLA). Almost all the Tour molecules have 
BLA together with donor and acceptor groups. In fact, 
the BLA which dimerizes the system have been shown 
to give rise to NDRM.. Interestingly, we find that donor 
(-l-e) and acceptor (— e) at some positions together with 
explicit dimerization also causes NDR and asymmetric 
I-V, very similar to the I-V shown in Fig. 3. The main 
point is that whether it is the explicit dimerization or 
two-sublattice structure, coupled with the voltage drop, 
it induces interchange of symmetry together with Landau 
quasi-degeneracy of the low-lying levels. 

To summarize, we find that the ratio i : £2 — ei is 



very crucial in determining the nature of the I-V char- 
acteristics. Three features, namely, the critical bias, the 
sharpness of the NDR peak and the extent of asymmetry 
in the I-V curves are sensitive to this ratio. We believe 
this to be the reason for experimentally observed NDR 
and asymmetry in Tour molecules containing NII2 and 
NO2 groups or NO2 group only, and its absence when 
there is only NII2 group or no substituents. Nitro group, 
having a very strong acceptor character tilts this ratio in 
favor of a sharp NDR peak and asymmetry, within the 
bias range considered. 

Although our model neglects some of the important 
issues in this field like electron-electron interactions. 
Coulomb blockade, and other many-body effects, it is 
nevertheless interesting to see how some experimentally 
observed features can be captured even by such sim- 
ple models. And a two-level model, although not well- 
supported by quantum chemical calculations, is not too 
far-fetched a description of the Tour molecules, in which 
the strength of the localized donor/acceptor reduces pro- 
gressively towards the chain ends. 

Our explanation for NDR encompasses the bias driven 
conformational change occurring in the system and the 
reduction of the acceptor to donor, where the roles of 
the donor and acceptor groups have been interchanged, 
as well as the bias driven changes in electrode-molecule 
coupling which would result in a change in the on-site 
energies of the components of the molecule. However, the 
common thread among all the known factors governing 
NDR requires further attention. 
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We also see an NDR peak at reverse bias region (at 
~ —2.61^ compared to ~ 1.6V in forward bias). This is 
because in the reverse bias, the charge movement is op- 
posed by the extra hopping term which arises due to the 
even number of sites in the chain. The difference between 
the critical bias in the positive and negative direction is 



governed dominantly by t and to a second order by the 
difference in the e. Note that, NDR on both sides of bias 
has been observed in an experiment with Tour molecule 
recently-. 



